Quantum and classical localisation, the spin quantum Hall effect and generalisations 
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We consider network models for localisation problems belonging to symmetry class C. This sym- 
metry class arises in a description of the dynamics of quasiparticles for disordered spin-singlet super- 
conductors which have a Bogoliubov - de Gennes Hamiltonian that is invariant under spin rotations 
but not under time-reversal. Our models include but also generalise the one studied previously in 
the context of the spin quantum Hall effect. For these systems we express the disorder- averaged 
conductance and density of states in terms of sums over certain classical random walks, which are 
self-avoiding and have attractive interactions. A transition between localised and extended phases 
of the quantum system maps in this way to a similar transition for the classical walks. In the case 
of the spin quantum Hall effect, the classical walks are the hulls of percolation clusters, and our 
approach provides an alternative derivation of a mapping first established by Gruzberg, Read and 
Ludwig, Phys. Rev. Lett. 82, 4254 (1999). 
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I. INTRODUCTION 

Localisation of a particle moving in a random environ- 
ment may occur both quantum-mechanically and with 
classical dynamics, but the phenomenon is very differ- 
ent in the two cases. In this paper we discuss a class 
of quantum-mechanical localisation problems for which 
some physical quantities can be expressed exactly in 
terms of averages taken in a classical counterpart. The 
equivalence holds despite the fact that interference effects 
dominate the behaviour of the quantum systems. 

Disordered quantum systems can in general be clas- 
sified according to their symmetries under time rever- 
sal and spin rotation. Three such symmetry classes are 
represented by the Wigner-Dyson random matrix ensem- 
bles, while an additional seven have been identified more 
recently. The models we study here belong to one of 
these additional classes, termed class C by Altland and 
Zimbauerfl One feature which distinguishes systems be- 
longing to each of the additional symmetry classes from 
those in the Wigner-Dyson classes is that they have a 
special energy in their spectrum, with eigenstates oc- 
curring in pairs either side of this energy. Some of the 
additional classes have realisations as Bogoliubov - de 
Gennes Hamiltonians for quasiparticles in disordered su- 
perconductors, where pairing interactions are treated at 
the mean-field level. Here, the special energy is the chem- 
ical potential in the superconductor and eigenstates are 
related in pairs by a particle-hole transformation. In par- 
ticular, class C arises for quasiparticles in a spin-singlet 
superconductor in which time-reversal symmetry is br 
ken for orbital motion but Zeeman splitting is negligible 
Since quasiparticle charge is not conserved in a supercon- 
ductor, experiments to investigate quasiparticle dynam- 
ics in these systems must probe thermal or spin trans- 
port. Moreover, as the characteristic features of the sym- 
metry class appear only close to the chemical potential, 



it is particularly gapless superconductors that are inter- 
esting: cuprate superconductors in the mixed state con- 
stitute a conspicuous example. 

The models we study are especially simple realisations 
of their symmetry class. They are obtained as gener- 
alisations of the network model originally introduced to 
describe localisation in the context of the integer quan- 
tum Hall plateau transition.! 2 ] Thus, they are formulated 
in the language of scattering theory and represent quan- 
tum particles, in general with TV-component wavefunc- 
tions, propagating on the directed links, or edges, of a 
lattice and scattering between links at nodes. The sym- 
metry of class C restricts N to even values, while our ap- 
proach requires that all nodes of the lattice have two in- 
going and two outcoming links. For these models, we are 
concerned with the density of states, obtained from the 
time-evolution operator, and with the disorder-averaged 
conductance of a finite sample, calculated from the Lan- 
dauer formula. In both cases, our starting point is an 
expansion for the Green function as a sum over Feynman 
paths. Our central result is that the terms in this sum 
which survive after disorder-averaging can be interpreted 
as self- avoiding classical random walks with attractive, 
short-range interactions. 

A particular network model from class C, in two i 
sions and with N = 2, has been studied previously.! 
It shows the so-called spin quantum Hall effect, having 
two insulating phases, with quantised values of Hall con- 
ductance differing by an integer, separated in the phase 
diagram by a delocalisation transition which is analogous 
to the quantum Hall plateau transition. In a remarkable 
paper, using supersymmetry to perform disorder aver- 
ages, Gruzberg, Read and Ludwigp (GRL) have shown 
that many physical quantities of interest in this model 
can be determined from the properties of the perimeters, 
or hulls, of classical percolation clusters in two dimen- 
sions. The approach we describe here provides an al- 
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ternative derivation of their results, using more elemen- 
tary, non-supersymmetric methods, as well as an exten- 
sion to other lattices, including ones in more than two 
dimensions and irregular lattices for which transfer ma- 
trix methods are inappropriate. It also extends to any 
even integer N. Our expressions give disorder-averaged 
physical quantities for the quantum system in terms of 
averages over classical random walks on the same lat- 
tice. In the case treated by GRL these walks are simply 
percolation hulls, for which many analytical results are 
available. By contrast, in the general case the proper- 
ties of the classical walks are not known. Nevertheless, 
the classical problem is much simpler than the original 
quantum problem, and we are able to construct further 
examples for which it is tractable. We remark that a dif- 
ferent type of connection between quantum Hall plateau 
transitions and percolation, based on the classical limit, 
has been discussed recently in Refs. || and |^. 

There are some important qualitative differences be- 
tween the properties of systems from Wigner-Dyson 
classes and those from the additional symmetry classes. 
In particular, while single-particle quantities such as the 
density of states are smooth functions of energy in the 
former case, in the latter case they may have singulari- 
ties at the special energy, which we take to be zero in the 
following. This is illustrated by previous results on be- 
haviour of mod^pfcom^la&s-jC, obtained using a variety 
of techniques .E300y'oE§E3 Random matrix ensembles 
with this symmetry, representing the zero-dimensional 
limit appropriate for quantum dots, have a density of 
states that vanishes quadratically in energy _aJ_ energies 
much smaller than the mean level spacing.EM Similar 
behaviour is expected for finite-dimensional systems if 
states are Anderson localised, on the grounds that ran- 
dom matrix theory should describe states within a local- 
isation volume.E3 Calculations for one-dimensio nal s ys- 
tems from class C, using either supersymmetryl 1 1 H 13 I or 
the DMPK equation,liil confirm this idea. 

Existing information on localisation in class C sys- 
tems is alsp, provided in part by renormalisation group 
treatments£3 of the appropriate non-linear sigma model 
at weak coupling, corresponding to weak disorder. These 
calculations identify two as the lower critical dimension. 
Thus, as in the Wigner-Dyson symmetry classes, it is 
only in more than two dimensions that a transition oc- 
curs between localised and metallic phases as a function 
of disorder strength, while in one dimension even weak 
disorder is sufficient to localise all states. In addition, 
for two-dimensional systems with broken time-reversal 
symmetry, including both ones from the Wigner-Dyson 
unitary class and ones from class C, a delocalisation tran- 
sition of the quantum Hall type is possible, and it is this 
transition that has-hewx-the focus of past work on class 
C network modelsU'DBtffl 

Many of these aspects, including the form of the den- 
sity of states in a localised phase and the possibility of a 
quantum Hall plateau transition, emerge naturally from 
the approach we describe here, which is presented as fol- 



lows. In Sec. II we introduce in detail the models that 



we are concerned with. In Sec. Ill we set out our general 
results, relating the density of states and average con- 
ductance for a network model to averages over certain 
classical random walks, and present proofs of these re- 
sults. We describe applications of these general results 
to the spin quantum Hall effect, to random matrix the- 
ory, and to localisation on a Cayley tree, in Sec. |iv| Open 
questions and future prospects are discussed in Sec.[v| 



II. MODELS 

We shall be concerned with models both for closed sys- 
tems and (in connection with the Landauer formula for 
conductance) for open systems, but we restrict definitions 
initially to closed systems. Consider a graph Q consist- 
ing of directed edges e connecting nodes n, each of degree 
four, with the restriction that at each node two directed 
edges enter and two leave. An ./V-component wavefunc- 
tion propagates on each edge. This propagation may 
be described by a unitary evolution operator U, which 
evolves the wave function one unit forwards in time, as 
the particle moves from a given edge to a neighbouring 
one. The evolution operator plays the same role in defin- 
ing the network model as does the Hamiltonian in the 
case, for example, of a tight-binding model. It has been 
discussed for the U(l) network model in Refs.[Tl] and [Ts|. 
In general, it is constructed from two ingredients. First, 
with each edge e is associated a unitary NxN matrix U e . 
This matrix specifies the (generalised) phase acquired on 
traversing the link. Second, with each node n is associ- 
ated an (S-matrix of the form 



S n = a 



cos 9 n sin 6 n 
— sin Or, cos 8„ 



(1) 



where 1 is the NxN unit matrix. This S'-matrix de- 
scribes scattering at the node from the incoming edges to 
the outgoing ones. If Q has E edges (and therefore E/2 
nodes), then U is an N x N matrix, with M = EN. It 
consists of E/2 blocks, each associated with a particular 
node and of size 2N x 2N. The block at the node n has 
the form 



U. 



1/2 

i 

U, 





1/2 



1/2 



u, 





1/2 



(2) 



where (1,2) and (3,4) label the edges which are respec- 
tively incoming and outgoing at this node. 

So far, the symmetry class of the network model has 
not been fixed, except that propagation along directed 
links breaks time-reversal symmetry n To identify net- 
work models from class C, one startsca from the defining 
pppperty of a Hamiltonian Ti with this symmetry, which 



TL* — —CTyTLdy , 



(3) 



where a y denotes the conventional Pauli matrix acting 
on spin variables and H* is the complex conjugate of 
H. Applying this to U, interpreted as U = e lH , the 
number of wavefunction components N must be even, 
so that the space of states on each link may be viewed 
as consisting of N/2 two-component subspaces, within 
which <jy operates. Then Eq. (0) becomes 



IA = (Ty U* <7y , 



(4) 



and from this an equivalent restriction follows on the 
edge phases, U e = o y U e *cr y , which are therefore unitary 
Sp(7V) matrices, equivalent for TV = 2 to SU(2) matrices. 

Randomness is introduced into these models via the 
edge phases. We take them to be independent random 
variables drawn from a distribution which is uniform on 
the invariant (Haar) measure of Sp(iV). The quenched 
average of a given quantity in the network model, denoted 
by (...), is the mean with respect to this measure. 

An open system is constructed from a closed system 
of this type by 'cutting open' a number of edges. The 
two halves of each edge cut into two in this way consti- 
tute one new edge directed into the system and one new 
edge directed out of the system. We may consider a con- 
ductance experiment between two 'contacts' by grouping 
a subset of these of the incoming edges {ei n } to form 
one contact and another subset of the outgoing edges 
{e ou t} to form the other. The transmission matrix t be- 
tween these two contacts is a rectangular matrix whose 
elements are (e ou t|(l — Z^) 1 1 ei rl ) . The spin conductance 
measured between the two contacts in units of (H/2) 2 /h 
is 
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TrtH . 



(5) 



from the multi-channel Landauer formula. 

Clearly, a great variety of specific models can be 
constructed within this framework, by making different 
choices for the graph Q and for the number of channels 
N. We defer discussion of particular examples to Sec. IV. 



III. GENERAL RESULTS 

In this section we state and prove our results for a 
general graph Q. We consider first the particular case 
of N = 2, and discuss the extension to general N in 



c 


s 


-s 


c 




FIG. 1: Decomposition of a given node. In the network 
model, S- matrix elements cos# n and ±sin# n are associated 
with the transitions (1,2) — > (3,4) as indicated on the left. 
Each decomposition of the node is then weighted with factors 
p n = cos 2 d n and 1 — p n = sin 2 8 n as indicated. 



Eq. (g) may be expanded as a sum over Feynman paths 
on Q which begin on e' and end on e: each path gives an 
ordered product of factors zllj along edges traversed by 
the path, weighted by appropriate factors of cos 9 n and 
±sin^ n for each node through which it passes. Alterna- 
tively, we may re- write Eq. (|6|) as 



G(e, e'; z) = -{e^U^l - ^V)" V) 



(7) 



obtaining instead a series convergent for \z\ > 1, involv- 
ing an ordered product of factors and an overall 
negative sign for each Feynman path. 

Our central result is an expression for the disorder- 
average, Sp(2) trace of the Green function, Tr (G(e, e, z)), 
in terms of classical paths. To state this result, we de- 
fine on the same graph Q a classical scattering problem 
as follows. Each node may be decomposed into two dis- 
connected pieces in two ways, viz. (13, 24) or (14, 23), as 
illustrated in Fig. [TJ. 

Theorem 1. The average Green function is given for 
\z\ < 1 by the generating function for the probability 
P(e; L) in the classical problem that the edge e belongs 
to a loop of a given length L. Explicitly 



Tr(G(e,e;z)} 



E 

L>0 



P(e;L)z 



2L 



For \z\ > 1 it is given instead by 

Tt(G{e,e;z)) = ^P{e;L)i 



-2L 



(8) 



(9) 



Sec. [HE 



L>0 



B. Density of states 



A. Green function, Feynman path expansion and 
classical walks 

The Green function for propagation from edge e' to 
edge e is (with N ~ 2) a 2 x 2 matrix 

G(e, e'; z) = (e|(l — zlA)~ Y \e') (6) 

where |e) is a state in the two-component space of wave- 
functions for a particle located on the edge e. For \z\ < 1 



The eigenvalues of the evolution operator U for a closed 
graph lie on the unit circle in the complex plane and may 
be written as exp(zej), with eigenphases — tt < ej < 7T 
for j = 1 . . . Af. These eigenphases are analogous for 
the network model to the energy eigenvalues of a system 
specified by its Hamiltonian. We define the density of 
states to be 

p(e) = ^3( £ - £j )). (10) 

j 
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A consequence of the symmetry of Eq. (||) is that p(e) 
Defining P(L), the edge-average of P(e; L), by 



P(L) = i^P(e;i), 



E 

e 

it follows from Eqns. (||) and (||) that 
1 



(11) 



2tt 



1 - ^ cos(2Le) 



L>0 



Conductance 



(12) 



The c lassical scattering problem introduced in 



Sec. Ill A may also be considered for an open system. For 



an open system with M cut links, each decomposition 
breaks Q into M directed paths which each run from an 
ingoing edge to an outgoing one, together with a number 
(possibly zero) of closed loops. Let P(e, e') be the prob- 
ability that a path runs from the ingoing edge e' to the 
outgoing edge e. 

Theorem 2. The disorder-average of the conductance 
defined in Eq. (1q) is given by 



(.9) =2 p ^ e 'y 



(13) 



eGl,e'G2 



where the sets 1 and 2 denote respectively: the edges 
incident on the first contact from Q; and those incident 
on Q from the second contact. 



D. Proofs 

It is useful to introduce for a closed system the resol- 
vent 



R(z) = ^TrG(e,e;z) 



(14) 



(where Tr again indicates an Sp(2) trace) and to gen- 
eralise this to the case where the parameter z takes in- 
dependent values z e on each edge e. The expansion of 
R({z}) as a sum over paths then yields a multinomial 
expression in all the z e . 

We now require two Lemmas about Sp(2) matrices. 
Lemma 1. If U G Sp(2), then its mean qih moment, 
(U q ), is zero unless q = or q = 2, in which case it takes 
the value 1 or — ^ 1 respectively. This can be shown using 
the representation U = exp(icm ■ a), where n is a unit 
real 3- vector and the as are the Pauli matrices, so that 



U q = cos qal + i sin qan ■ a . 



(15) 



The result now follows from the observation that the in- 
variant measure for Sp(2) is the uniform measure on the 



group manifold S3, and therefore has the form J dU — 
7r ~ 1 IqQ ~ cos2a)da J dri. 

Lemma 2. If G is a real linear combination of Sp(2) 
matrices, it is itself proportional to an Sp(2) matrix, with 
a real scalar constant of proportionality. This follows 
directly from the above representation for each matrix. 

The main argument in the proof of theorem 1 now pro- 
ceeds as follows. As discussed above, for each realisation 
of the randomness, R({z}) is a sum over closed directed 
paths on Q . For a particular path, each link or node may 
be traversed an arbitrary number of times. We first show 
that it is sufficient, in calculating the mean (R({z})), to 
restrict this sum to those paths which traverse each edge 
exactly twice or not at all. Let us consider the sum of all 
paths in which a particular edge e is traversed exactly q 
times. This has the form 



zfTr [U e A(e, e) 



(16) 



where A(e, e) denotes the sum over all weighted paths 
which begin and end on e, but do not themselves traverse 
e. By Lemma 2, this is proportional to an Sp(2) matrix,E£l 
so it may be written as 



A(e,e) = \A(e,e)\A(e,e) 



(17) 



where \A(e, e)| is real and A{e, e) £ Sp(2). Defining U' e 
U e A(e,e), (111) may e written 



zlTrU?\A(e,e 



(18) 



The invariant integration over U e is equivalent to that 
over U' e , so that, by Lemma 1, the result will vanish un- 
less q = or q — 2. Since (R({z})) is a multinomial 
expression in the parameters {z}, the argument may be 
applied to each edge in turn to show that the only al- 
lowed powers of any z e entering this expression are or 
2. This establishes the first part of the proof. 

For paths in which each node is visited only or 2 
times, the main result follows immediately. For such a 
path must traverse a closed loop in the decomposition 
of G exactly twice. At each node there will be factors of 
cos 2 8 n or (± sin6*„) 2 , giving precisely the correct weight- 
ing for this loop to appear in the decomposition. The 
product of Sp(2) matrices along the loop will have the 
form 



Tr({U 1 U 2 ...)(U 1 U 2 ...)) 



(19) 



Thus, defining U[ = U1U2 . ■ ■ , this is equivalent to av- 
eraging U[ , which gives — ^ by Lemma 1. Finally, the 
trace gives a factor of 2. 

Further effort is required to treat paths which visit 
some nodes more than twice. A little thought shows that 
such a node must be visited exactly four times, enter- 
ing and leaving exactly twice along each directed edge, 
if the contribution to R(z) is not to vanish on averaging. 
Consider such a node n, and label the incoming and out- 
going edges as (1, 2), (3,4) respectively (see Fig. Q). We 
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show that the sum over all paths visiting this node four 
times may be written in terms of the two ways of decom- 
posing this node, with precisely the correct weights. Let 
A(i,j) be the sum over all paths from edge j £ (3,4) to 
edge i 6 (1, 2), which do not use any of these four edges. 
The sum over paths visiting the node four times may 
be decomposed into eighteen different contributions, de- 

I 



and 



where we have introduced the shorthand s = siriO n and 
c = cos 9 n . (Note that ( pi] ) and ( |22| ) gi ve separate con- 
tributions to R({z}); similarly for ( p4|) and (p5|)). As 
before, write A(i, j) = \A(i,j)\A(i,J) and note that, by 
a change of integration variable, (GO) is, on averaging, 
equivalent to the average of U 2 (that is, —5), multiplied 
by \A(1, 3)| 2 | J 4(2,4)| 2 s 4 . By a similar argument, each of 
( f2l| ) and ( |22] ) are equal to the average of U 2 U 2 (that 
is, (-|) 2 ), multiplied by I A(l,3)| 2 |A(2,4)| 2 c 2 (-s 2 ). The 
total of the (§(]) - © is therefore 

- L4(l, 3)| 2 L4(2,4)| 2 ( s 4 + 2(-I)c 2 (- s 2 )) 
= -|A(l,3)| 2 |A(2,4)| 2 sin 2 0„ 

I 



In addition, there are six equal contributions in which 
the factors are cyclically permuted so that each be- 
gins Uj.A(l, 4) . . . . Each term is proportional to 
|A(1,3)||A(2,4)||A(1,4)||A(2,3)|. The remainder of the 
expressions may be simplified, for example, by redefin- 
ing ?7iA(l,3) = U{ and U 2 A(2,4) = U' 2 . This reduces 



pending on the order in which the edges are visited. Each 
contains a product of four factors A(i, j) as well as Sp(2) 
matrices Ui and Uj . There are in fact two types of contri- 
bution: in the first type, two of the A(i,j) appear twice 
and the others not at all; while in the second type all 
four A(i,j) appear once each. There are six of the first 
type and they have the form 



(20) 
(21) 
(22) 



(23) 
(24) 
(25) 

I 

Similarly, (p3]-p5[) sum up to 



— 1-4(1, 4) | 2 \A(2, 3) | 2 cos 2 e n (27) 



Now consider the other twelve contributions, in which 
A(l,3), A(2,4), A(l,4) and A(2,3) each appear exactly 
once. Six of these are 



(28) 
(29) 
(30) 

(31) 
(32) 

(33) 

I 

(IS to 

Tr U[U 3 U{BU 4 U 2 U 4 U^CU 3 c 4 (34) 
Tr U[U 3 U 2 CU 3 U 2 U 4 U[BU 4 s 4 (35) 
Tr U' 1 U 3 U[BU 4 U 2 CU 3 U 2 'U i c 2 (-s 2 ) (36) 
Tr U[U 3 U' 2 C\J 3 U[BlJ 4 U l 2 U 4 c 2 {~s 2 ) (37) 
Tr U[U 3 U 2 U 4 U[BU 4 U 2 CU 3 c 2 (-s 2 ) (38) 
Tr U{U 3 U 2 U 4 U 2 CU 3 U[BU 4 c 2 {-s 2 ) (39) 



Tr C7iA(l, 3)U 3 U 2 A{2, 4)U i U 1 A(l,3)U 3 U 2 A(2, 4)C/ 4 s 4 

Tr U X A(1, 3)U 3 U x A(l,3)U 3 U 2 A(2, 4)U 4 U 2 A(2, 4)U 4 c 2 (-s 2 ) 

Tr U X A(1, 3)U 3 U 2 A(2, 4)U 4 U 2 A(2, 4)U 4 U X A(1, 3)Z7 3 c 2 (-s 2 ) 



Tr U X A(1, 4)U 4 U 2 A(2, 3)U 3 U l A{l,4)U 4 U 2 A(2, 3)U 3 c 4 
Tr U X A(1, 4:)U 4 UiA(l,4)U 4 U 2 A(2, 3)U 3 U 2 A(2, 3)U 3 c 2 {-s 2 ) 
Tr U X A(1, 4)U 4 U 2 A(2, 3)U 3 U 2 A(2, 3)U 3 U X A(1, A)U 4 c 2 (-s 2 ) 



Tr U X A(1, 3)U 3 U x A(l,4:)U 4 U 2 A(2, 4)U 4 U 2 A{2, 3)U 3 c 4 
Tr U X A(1, 3)U 3 U 2 A(2, 3)U 3 U 2 A(2, i)U 4 U x A(l, 4)C/ 4 s 4 
Tr U X A(1, 3)U 3 U X A(1,4)U 4 U 2 A(2, 3)U 3 U 2 A{2, 4)U 4 c 2 {-s 2 ) 
Tr U X A(1, 3)U 3 U 2 A(2, 3)U 3 U X A(1,A)U 4 U 2 A{2, 4)U 4 c 2 (-,s 2 ) 
Tr U X A(1, 3)U 3 U 2 A(2, 4)U 4 U X A(1,A)U 4 U 2 A{2, 3)U 3 c 2 (-,s 2 ) 
Tr U X A(1, 3)U 3 U 2 A(2, 4)U 4 U 2 A(2, 3)U 3 U X A(1, 4)U 4 c 2 {-s 2 ) 



where B = A(l, 3)- 1 A(l, 4) and C = A(2, 4)- 1 A(2, 3). 
By making suitable changes of integration variables, 
as before, we find that (|34|) and ( |35| ) give factors of 
(-i) 2 Tr (BC), while © - (£§ give (-±) 3 Tr (BC). The 
sum of all twelve such contributions is therefore 

1. 



X = -Tr(BC)\A(l,3)\\A(2,i)\\A(l,A)\\A(2,3)\ 
x (c 4 + .s 4 + 2c 2 s 2 ) 
= iTrA(l,3)U(l,4)A(2,4)U(2,3) 



(40) 



The important feature of this result is that it is indepen- 
dent of 9 n . It may be written, trivially, as 



X = Pn X + (1 - Pn )X 



(41) 



while the first six contributions have the form (as given 
in H and |27|) 

-p n L4(l,4)| 2 L4(2,3)| 2 - (1 -p„)L4(l,3)| 2 L4(2,4)| 2 

(42) 

Therefore we can obtain the same total result, after av- 
eraging over U\, . . . , t/4, if we decompose the node n in 
each of the two possible ways, and weight the two de- 
compositions with probabilities p n and 1 — p n respec- 
tively (see Fig. nl). For in the first case, (13,24), we 
find -|A(1,4)| 2 |A(2,3)| 2 +X, while in the second case of 
(14, 23), we get -\A(l, 3)| 2 |A(2, 4)| 2 + X. 

Now we may simply go through Q, decomposing it 
node by node. The result is that (R({z})} is given by 
the weighted sum over the same quantity calculated on 
graphs arising from all possible decompositions of Q\ such 
graphs consist of closed loops with no remaining nodes. 
In each decomposition, a given edge e belongs to just one 
loop, and, for that loop, of length L say, the contribution 
to Tr G(e, e; z) is just —z 2L . Thus the mean Tr (G(e, e; z) 
is given by the average of this quantity over the ensemble 
of loops corresponding to the decomposition of Q. This 
completes the proof of Theorem 1. 

So far we have considered unitary evolution on a closed 
system. However, since the proof is constructed to work 
for arbitrary values of the parameters {z e }, it generalises 
straightforwardly to an open system where probability 
is not conserved at the nodes which are connected to 
external leads. For example, if in Fig. |l|, the edges (1, 3) 
correspond to external leads which carry no current, then 
we should consider only those Feynman paths through n 
which pass directly from 4 — > 2, with amplitude cos# n . 
In the above argument, this may be taken into account by 
regarding this node as a single edge, carrying an Sp(2) 
matrix U4U2, and with fugacity ziz^cosQ n . With this 
modification, the proof of Theorem 1 goes through as 
before. 

Now consider conductance measurements on a open 
graph Q, with external contacts {e ln } and {e ou t}, as de- 
scribed in Sec. II. Consider the graph Q' formed from 
Q by joining particular outgoing and incoming external 




s 



DC 



FIG. 2: Implementation of a link in the Sp(4) model. It con- 
sists of pairs of incoming and outgoing Sp(2) links, connected 
through a random S-matrix which mixes the two channels. 
After averaging, this is equivalent to the decomposition shown 
on the right hand side, equally weighted. 



ed ges e ou t and e- ln to create a new internal edge e. Ob- 
serve that, in calculating (G(e, e; z)), the edge e is always 
traversed exactly twice. Therefore 



(G{e ou t,e in ;z) 2 )g = (G(e,e;z))c 



(43) 



Since G(e ou t,ei n ) is given by a sum over Sp(2) ma- 
trices, by Lemmas 1 and 2, (G(e ou t, e; n ; z) 2 ) = 
— |(|G(e ou t, e; n )| 2 ). Now apply Theorem 1: it follows 
that (|G(e ou t, ei n )| 2 ) is twice the probability that ei n and 
e ou t are connected by a path in the decomposition of Q. 
Summing over all the edges e- m and e ou t as required by 
the Landauer formula (|5|) then gives the result stated in 
Theorem 2. 



E. General N 

Our approach to the Sp(7V) model with even N > 2 
is based on the fact that a general Sp(iV) rotation can 
be written as a product of rotations derived from a lim- 
ited number of generators. Moreover, under rather weak 
conditions, the probability distribution of such a product 
containing random factors will converge as the number 
of factors increases to the Haar measure on Sp(N): the 
requirement is simply that no subspace is left invariant 
by the ensemble of rotations. 

We therefore build the Sp(7V) model on Q by taking 
N/2 copies of the Sp(2) model, each individually defined 
on Q, and coupling the copies together. This coupling 
takes the form of a product of many non-commuting 
matrices. Successive factors in the product are of two 
kinds. One kind is block-diagonal with N/2 random 
Sp(2) blocks, resulting in intra-copy scattering; the other 
kind must produce inter-copy scattering and may be cho- 
sen non-random. It is convenient to restrict the inter- 
copy scattering to be between corresponding links in each 
copy. For example, in the case of Sp(4) each pair of cor- 
responding links is coupled as shown in Fig. |2|. The evo- 
lution matrix for this pair of links consists of a product 
S of a large number of factors, each of the form 



J_ [Ui 
V2 \ U 2 



1 1 

-1 1 



(44) 



where Ui and U2 are random Sp(2) matrices, chosen in- 
dependently for each term in the matrix product. 
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The advantage of this choice is that we may immedi- 
ately apply our main theorem to show that, after averag- 
ing over the random Sp(2) matrices, the mean density of 
states and mean conductance are given in terms of a de- 
composition of each node as before, and a decomposition 
of each link as shown in Fig.||, with equal probabilities 
for each term. This construction may be generalised to 
Sp(iV) for arbitrary even TV: the result is that, in the 
classical system, the N/2 exit trajectories are a random 
permutation of the incident ones, with all permutations 
having equal weight, independently for each edge. 



IV. APPLICATIONS OF GENERAL RESULTS 

There are three obvious classes of behaviour possible 
in the classical scattering problem which we have arrived 
at, corresponding to loops, or closed classical walks, that 
are localised, extended or critical. 

By localised classical walks, we mean that only a van- 
ishing fraction of walks have infinite length, and that 
the number of walks longer than a characteristic size £ 
decreases rapidly with £. One might have, for example, 
P(L) <~ cxp(— £/£) for L ^> £. Then the fraction of loops 
longer than L, given by 



/(L) = l-£p(Z), 



(45) 



1=1 



approaches zero as L — > oo. As a consequence, p(e) — > 
as e — > 0. Introducing the mean square length of loops, 



(L 2 ) = £ P(L)L 2 , 



(46) 



L>0 



then provided (L 2 ) is finite, p(e) vanishes quadratically, 
varying as 



(47) 



This is the behaviour expected in the localised phase 
of the quantum problem. Moreover, if classical walks 
have a characteristic size £ and if Q is embedded in Eu- 
clidean space, then we expect the conductance to de- 
crease rapidly with increasing contact separation, for sep- 
arations larger than £. 

Alternatively, classical walks are extended if a finite 
fraction have infinite length, so that Xvca.^^^ f(L) > 0. 
In this case, p(0) — (2n)^ 1 f(oo) is non-zero, which we 
expect in the extended phase of the quantum problem. 
Clearly, the form of p(e) at small e is determined by that 
of f(L) at large L: if f(L) - /(oo) - L~ x with < x < 2 
then p(e) — p(0) oc \e\ x , while if f(L) — /(oo) falls faster 
than L~ 2 at large L, then p(e)_— p(0) ~ e 2 . It is an in- 
teresting consequence of Eq. ( |12j ) that p(e) is bounded, 
and so a divergence in the density of states, as occurs 
for example in class D localisation problems, is excluded 
in class C network models, whatever the choice for Q. 
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FIG. 3: The network model for the spin quantum Hall effect, 
defined using the L-lattice. Scattering amplitudes of ± COS C7 n 
and ±sin# n are indicated by ±c and ±s. 



In the extended phase for a system in d-dimensional Eu- 
clidean space, Ohm's law dictates that the two-terminal 
conductance should vary with the separation, Zy, between 

terminals and their cross-sectional area, l±, as Ij^ 1 /l\\ for 
large Zy, l±; this places strong constraints on behaviour 
in the corresponding ensemble of classical walks, suggest- 
ing that, on large distance scales, they should behave like 
free random walks. 

Classical walks that are critical have no characteris- 
tic loop size, and a vanishing fraction that are of infinite 
length. Then /(oo) = and a possible behaviour for 
P(L) is P(L) - L~y at large L, with y > 1. The result- 
ing quantum density of states has the critical behaviour 
p(e) ~ |e| a_1 for small e. The conductance of such a 
system is expected to depend on the geometry of sample 
and contacts, but to be unchanged under a rescaling of 
all spatial dimensions. 

The task that remains, given a particular model for 
quantum localisation, specified by the graph Q and values 
for the node probabilities p„, is to determine behaviour 
of the corresponding classical walks. In general, this re- 
mains a challenging open problem with connections to 
previously-studied random- walk problems which we sum- 
marise in Sec. |v|. In particular instances, however, rele- 
vant properties of the classical walks can be calculated; 
we describe three examples in the remainder of this sec- 
tion. 



A. Spin quantum Hall effect 

A two-dimensional model exhibiting the spin quantum 
Hall effect is obtained by taking Q to be the L-lattice, 
illustrated in Fig. |[ As GRL have shown, the two possi- 
ble classical decompositions of a node may be associated 
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with the presence or absence of a bond, with probabilities 
p and 1 — p, between neighbouring sites on an associated 
square lattice. This associated lattice is rotated by 45° 
relative to the L-lattice, and has a larger lattice spacing 
by a factor \[2. In this way, closed loops of the classical 
problem form interior or exterior hulls of bond percola- 
tion clusters on the larger lattice. It is known that such 
loops are finite with characteristic size £ except at the 
critical point, p = p c , which for bond percolation on the 
square lattice occurs at p c = 1/2. On approaching the 
critical point, £ diverges as £ ~ \p — p c \~ v with v = 4/3, 
while at the critical point the distribution of hull lengths 
is P(L) ~ L~ v at large L, with y = 8/7. In this way 
one finds that the quantum localisation length diverges 
with the same exponent value, v = 4/3, as the plateau 
transition is approached, and that the density of states 
p(e) varies for small e as p[e) ~ e 2 in the localised phase, 
and as p(e) ~ H 1 / 7 at the critical point. 



o 



B. Random matrix theory 



The simplest application of our discussion of Sp(iV) 
models with general N is to random matrix theory. To 
this end, we take Q to consist of a single edge closed on 
itself. Then the evolution operator U is a random Sp(iV) 
matrix, chosen with the Haar measure. In this case the 
density of states p(e) is the eigenphase density for the 
Sp(iV) random matrix ensemble- which has been deter- 
mined previously by Zirnbauer£9 using supersymmetry 
methods. 

Applying to this problem the approach we have de- 
scribed in Sec. [HE , we must consider N/2 copies of an 
edge, which are closed by connecting outgoing ends to 
a permutation of ingoing ends. For example, the case 
N = 4 corresponds to joining opposite ends of the two 
possibilities shown in Fig. |J thus producing with equal 
probability either two loops, each of length 1, or a sin- 
gle loop of length 2. For general N, all possible lengths 
L up to N/2 are possible, with equal probability. Thus 
P(L) = 2JN for 1 < L < N/2, and P(L) = otherwise. 
Using (|l2|) we thence find 



N + l 



2ttN 

in agreement with Rcf. nd. 



1 



sin(iV + l)e 
(N + l)sine 



(48) 



C. Cay ley tree 

A solvable model which illustrates each of the types of 
behaviour - localised, extended and critical - for classical 
loops is based on the geometry of the Cayley tree. Specif- 
ically, we take Q to be a graph of the type illustrated in 
Fig 4; the U(l) network model on such a tree has been 
studied previously in Ref. ^TJ. We restrict attention to 
N = 2 and consider first coordination number q — 3. As 




FIG. 4: The graph Q for network model based on a rooted 
Cayley tree. An example with coordination number q = 3 
and four generations is shown; system size is increased by 
increasing the number of generations. S-matrix amplitudes at 
nodes are cos 8 to remain at the same generation, and ± sin 8 
to change generations. 



previously, we define P(e; L) to be the probability that 
the edge e of the root, far from the surface of the tree, 
lies on a given closed loop of length L. 
Then 



P(e; 3L) = P 2 6 LA + 2p{l - p)P(e; 3L - 3) 



L-2 



-(1 - pf p {e; 3m)P(e; 3L - 3m - 3) . 



(49) 



This can be solved by using the generating function 
G(z) = J2 r >i zr P( e > 3r): we obtain 

G(z) = zp 2 + 2zp{\ - p)G(z) + (1 - pfzG{z) 2 . (50) 



We take the branch of G(z) such that G(0) = 0; thus 
1 



G(z) = 



2z{l-pf 



(l - 2zp(\ -p)- v/l-4zp(l-p)) 

(51) 



When z — 1 , yl — 4zp(l — p) is equal to 1 — 2p for p < 
l/2and2p-lforp> 1/2. Therefore for p > 1/2, G(l) = 
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1, so that all walks are localised. The characteristic loop 
size is given by 



< L > = 3 



dG 
dz 



2p- 1 



(52) 



However, when p < 1/2, G(l) = [p/(l ~ p)} 2 < 1 . Thus 
in this case a proportion 1 — G(l) of the walks do not 
close; they are extended. For such walks, (^9|) is invalid 
since it describes only closed walks of a finite length. For 
p = 1/2 and large L 6 Z, we find the critical behaviour: 
P(e; 3L) ~ . The corresponding density of states in 
the quantum system (defined locally on an edge far from 
the surface, to avoid surface effects) has the behaviour 



p(e)oc^|ep p 

constant p < 



P > 2 

_L 

2 
1 

2 



(53) 



We can extend this analysis to consider coordination 
number q > 3; then 



G(z) = z[p+(l-p)G(z)} 



9-1 



(54) 



Thus z b = (1 - - i) 1-9 !^ - 2)/p] 9 ~ 2 locates the 

unique branch point of G(z). The complex roots of the 
two branches of G(z) corresponding to G(l) = 1 (all 
paths of finite length) and G(l) = W (where W is the 
fraction of paths which are infinite) , coincide at p = (q — 
2)/{q - 1). We conclude that for p > (q - 2)/{q - 1) all 
walks are localised, whereas for p < (q — 2)/(q — 1) there 
exists a fraction W of extended walks. In addition, we 
find that the critical behaviour P(e; qL) ~ holds for 
all q > 3. In the localised phase, the typical loop length 
is given by 



< L >-- 



p-(l-p)(q-2) 



(55) 



Thus, with critical value p c = 1/2 replaced by p c = 
(q — 2)/{q — 1), the behaviour of the density of states 
is the same as in the q = 3 case. The critical exponents 
are therefore universal and independent of q. In fact, we 
note that there is a one-to-one correspondence between 
walks on the Cayley cactus and percolation clusters on 
the Cayley tree; therefore the critical exponents are re- 
lated to those of percolation. 



V. DISCUSSION 

In this section we outline directions for future work. 
These are of two kinds. First, one can imagine a vari- 
ety of models which are likely to exhibit phenomena not 
shown in the examples treated in Sec. |v| Second, one 
can hope to make use of connections between the classical 
random walks that arise in our approach and statistical 
problems studied previously. Further work on both these 



aspects is in progress.0El In general, the equivalent clas- 
sical problems correspond to self-attractive random walks 
of various kinds. They are attractive because the weight 
for passing through a given node twice is p, rather than 
p 2 . However, the actual behaviour is expected to de- 
pend very much on the individual lattice. Many of these 
problems correspond to the classical scattering of light by 
random arrays of mirrors: for example, in each decom- 
position of Fig. [j], a two-sided mirror can be placed so 
that the classical trajectory reflects off it. Such problems 
have been studied extensivelyc3o'E§E3. They may also 
be realised in terms of history- dependent kinetic random 
walks, in the sense that the walker may be thought of 
as placing a mirror at random the first time it reaches 
a node: if it revisits the node, however, the mirror is 
already in place. This puts such models in the class of 
so-called 'true' self-avoiding walks, which have been stud- 
ied using field-theoretic RG methods in the limit of weak 
scatteringo. Interestingly, such studies indicate a critical 
dimension of two, just as for Anderson localisation. 

A two-dimensional model with behaviour which we ex- 
pect to contrast with that found in the spin quantum 
Hall effect can be obtained by taking Q to be the Man- 
hattan lattice, illustrated in Fig. ||, in place of the L- 
lattice. The crucial distinction is that, for the network 
model on the L-lattice, two distinct phases can be identi- 
fied from the different nature of edge states in the limits 
p — > and p — > 1, while for the Manhattan lattice there 
are no edge states in either limit. Equivalently, the Hall 
conductance of the model on the L-lattice may be non- 
zero, being determined on short distances by the value 
of p, but quantised at large distances, while on the Man- 
hattan lattice the Hall conductance always has average 
value zero. Consequently, one expectSrijom the scaling 
flow diagram for quantum Hall systemsa and the renor- 
malisation group calculations ior the class C non-linear 
sigma model at weak couplingt3 that states of the Man- 
hattan model should be localised for all p > 0. Since 
classical trajectories on the Manhattan lattice cross at 
nodes with probability p, they are not the hulls of perco- 
lation clusters. Nevertheless, it is straightforward to use 
bond percolation on a lattice which has a lattice spacing 
larger by a factor \[2 and is rotated by 45° relative to the 
Manhattan lattice to set an upper bound on loop sizes: 
this is sufficient to prove localisation for p > 1/2, and for 
p = 1 classical loops are simply the elementary plaque- 
ttes of Q. Conversely, at p — classical trajectories are 
simply straight lines, and so one expects the localisation 
length to diverge as p — > 0. 

A second variant on the spin quantum Hall effect can 
be constructed by considering the Sp(iV) model on the 
L-lattice, but with N > 2. It is natural to anticipate 
from the_scaling flow diagram proposed for quantum Hall 



systemsES that this model should have N/2 delocalisation 
transitions as the node parameter p increases from p — 
to p = 1. Between transitions, states are localised 
and the Hall conductance is quantised. As a transition 
is approached, the localisation length diverges, and on 



10 



FIG. 5: The Manhattan lattice. In the network model, S- 
matrix elements of cos 9 and ± sin 6 are associated respectively 
with 90° turns at nodes and with propagation in a straight 
line. 



passing through a transition, the number of edge states 
and the quantised value of the Hall conductance both 
change by two. It remains a challenge to understand 
how or whether such behaviour arises in the language 
of classical walks, and might be interesting to relate the 



large- TV limit of the lattice problem to the field theory 
discussed in Ref. |l6| 

A third direction is to consider models defined on 
graphs in three or more dimensions. A three-dimensional 
version of the U(l) network irraeLel has been studied pre- 
viously using a layered systemE2l and models with cubic 
symmetry may also be constructed. For systems in three 
and higher dimensions one expects that both localised 
and metallic phases should be possible, each existing over 
a range of values for node parameter, p, with a transition 
between the two at a critical value p — p c . Since in two 
dimensions it is known that the theta-point transition 
between collapsed and swollen phases of a self-attracting 
polymer chaip-,has the same exponents as those of per- 
colation hullsE3, it is tempting to suggest that in higher 
dimensions this delocalisation transition might also be in 
the theta-point universality class. If so, we would ex- 
pect mean-field behaviour, with logarithmic corrections 
for d = 3. 
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